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For N >5 there is a first order bulk transition that cleanly separates the strong and weak coupling 
regimes of SU(A0 lattice gauge theories with the plaquette action. We find that in this case the 
calculated string tension can be readily fitted throughout the weak coupling region by a standard 3- 
loop perturbative expression modified by lattice spacing corrections of the expected form. While 
our fits demand the presence of the latter, they are not constraining enough to tell us which of the 
various bare coupling schemes is a 'good' one, in the sense that terms in the j3 -function beyond 
3-loops are indeed negligible (in the relevant range of scales). To resolve this ambiguity we work 
in SU(3), using the Schrodinger Functional coupling scheme as a benchmark, and find that the 
Parisi mean-field improved coupling scheme matches it very well. Using the latter scheme, we 
have fitted the values of the string tension a 2 a that have been calculated for 2 > N > 8, to obtain 
A Ms/\/ a = 0.503 (2) (40) +0.33 (3) (3)/A^ 2 for N > 3, where the first error is statistical and the 
second is our estimate of the systematic error from all sources. 
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1. Introduction 

Consider SU(/V) lattice gauge theories with the standard plaquette action: 

Z= |n^expj-/3£ {l-^ReTr^, jj (1.1) 

where U p is the ordered product of the SU(A0 matrices around the boundary of the plaquette p. 
The parameter j3 is the inverse bare coupling, and this defines a running coupling on the scale a in 
what one can call the 'Lattice' coupling scheme: 

2N 

(1.2) 



It would be convenient to be able to determine a in units of a physical quantity, say the string 
tension a, from the value of g\(a) using a weak coupling expansion of the form: 



a^{a) ~ (l+ca 2 o + 0(a 4 ))F PT (g 2 L (a)) (1.3) 

where FpT{g\(a)) is obtained by integrating the continuum j8 -function at some (practical) order in 
perturbation theory. The additional factor containing an 0(a 2 ) correction with coefficient c ~ 0(1) 
must be there ^ since if we were to use some other physical quantity \l ! in place of \l = \Ja we 
would in general have 

with c' ~ 0(1), not to mention any 0(a 2 ) corrections from the j3 -function on the lattice. 
There are two well-known problems with implementing this: 

• g\ is a poor expansion parameter, as indicated by 

which implies that the L scheme will have large higher order terms in the /3 -function (assuming 
that the MS scheme is a 'good' one and does not); 

• it is not clear at what j8 we should expect such a weak coupling expansion to begin to work well, 
since SU(3) has a smooth strong-to-weak coupling crossover where 

powers in j8 — > powers in — , (1.6) 

and this makes it hard to evaluate the relative merit of an 'improvement' to the lattice-scheme from 
an apparent success in fitting a wider range of bare couplings. 

In this talk we describe the following strategy to resolve these two obstacles. First we use 
the fact that for SU(/V > 5) there is a first order 'bulk' transition [0|, that separates the weak and 
strong coupling ranges, thus removing the ambiguity of where one might expect a weak coupling 
expansion to be applicable. (Just like the Gross-Witten transition [Q] in D = 2.) While this en- 
ables us to quantify the importance of retaining 0(a 2 ) lattice corrections, it does not enable us to 
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usefully discriminate between various bare coupling schemes which lead to quite different values 
for Ajfg/y/c Presumably some have large higher order corrections in their /3-function and so are 
'bad'. To determine which of the schemes are 'good' ones we return to SU(3) and make use of 
the accurate calculation of the running coupling in the 'Schrodinger functional' (SF) scheme, that 
covers an energy range comparable to that of experiment, i.e. up to ~ Mz, and with appreciably 
smaller errors [|J]. We shall use this scheme to obtain, from the values of a/rQ calculated in [||] 
the continuum value of r^AsF and hence of foA^. We compare this to what one obtains with var- 
ious improved bare coupling extrapolations, and find that the original Parisi mean-field improved 
scheme [|]] closely matches the SF result. We simultaneously perform a comparison with the SF 
scheme that does not involve the calculation of a physical quantity and therefore can be carried out 
to much weaker coupling. This also points to the 'goodness' of the mean-field scheme. Motivated 
by this we use the latter scheme for N / 3 to obtain continuum values for Aj^/ ^fa for all N, and 
in particular for N — > °°. 

In this talk we present a brief summary of our work: details, including estimates of the various 
systematic errors, will be published elsewhere 



2. Lessons from larger 

In Fig. 1 we see the bulk transition, and its large metastability region, for SU(8). 
1.2 

ay/a 
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Figure 1: The SU(8) string tension versus the inverse lattice coupling, including the region of the first order 
'bulk' transition between strong and weak coupling. Values o are obtained coming from strong coupling, 
while the values • are obtained coming from weak coupling. 

In Fig.2 we show a fit to the weak coupling branch, all the way to the extreme metastability 
edge, using 
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where the scheme being used is the Parisi Mean Field Improved coupling 
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(2.1) 



(2.2) 



where U p is the plaquette variable. In eqn(2.1) the terms that involve only j3o and j8i constitute the 
exact 2-loop continuum result. (That is to say, it is the exact result when j3 ; >2 = 0.) We present 
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gj(a)N 




Figure 2: The 't Hooft coupling, defined from the mean-field improved lattice bare coupling as a function 
of the scale a in SU(8). Shown is the 3-loop perturbative running modified by a 0(a 2 ) lattice correction. 

the 3-loop contribution as a power series in g 2 . We note that although the coefficient c is actually a 
power series in gj, within our accuracy it suffices to treat it as a constant. 

The fit to SU(8) has c = 1.1 8 ±0.04 confirming the need for 0{a 2 ) corrections with coefficients 
of 0(1). However if we vary the perturbative coupling scheme we find that the range and accuracy 
of our calculations does not discriminate usefully between them. 

Comparing the values of gj(a)N for various N at fixed a^Jo, shows good evidence for a large- 
N j6 -function with very small corrections except at coarse lattice spacings. Thus it makes sense to 
take what we learn in SU(8) as a basis for treating other N, in particular SU(3). Performing fits with 



eqn(2.1) in SU(3) one sees in Fig || that these are only acceptable for j3 > 5.9, i.e. a^/o < 0.25, 
in contrast to the range a^Jo < 0.42 for SU(8). For SU(2) the range is even more limited, i.e. 
a^/o < 0.18. This shows explicitly how the smoothening of the strong to weak coupling transition 
means that one has to go to much smaller values of a to be able to use weak coupling expansions. 



gj(a)N 




ay/a 



Figure 3: The 't Hooft coupling, defined from the mean-field improved lattice bare coupling as a function 
of the scale a in SU(3). Shown is the 3-loop perturbative running modified by a 0(a 2 ) lattice correction. 



3. Choosing a good coupling scheme 

To choose a good bare coupling scheme s, we calculate A s / /I and hence Ajj^/jj., within various 
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such schemes (for some physical mass }i) and find which scheme produces values that agree with 
what we obtain using a 'reliable' lattice coupling scheme. For the latter we take the Schrodinger 
functional scheme of the Alpha Collaboration which for SU(3) [Q] covers a range of energy scales 
comparable to that covered by experimental measurements, and does so with greater precision. 
(Compare Fig.4 of [Q] with Fig.10 of The coupling g 2 SF has been calculated for a wide variety 
of values of j8 on scales la(p) where typically / = 6 to 12. We then take the calculated values of r^/a 
in ^ and interpolate these to the values of j8 at which g 2 SF (la) has been calculated. (Interpolating, 
unlike extrapolating, is a well controlled process.) We then fit using 

la 1 ( . ? r a 2 1 

'o 



1 + c" p Id 

r (a) r A SF V r ^ r W 



1 \W -fa&M 



\Po Pogj F {la)J 

Here there are two lattice spacing corrections. The usual 0(a 2 ) term arises from corrections to 
ro(a) etc. while the 0(l/l p ) term arises from lattice corrections to g 2 SF {la) on the scale / x a. We 
perform fits with both p = 1 and p = 2 taking the difference as part of our estimate of the systematic 
error. We obtain 

L_ = 3.2(1) — > ro A m = 0.640(20) (3.2) 
roJ^sF 

We now repeat this calculation using several lattice bare coupling schemes in fits of the form 
in eqn(2T) but with a^a(a) replaced by a/ro(a). For the Parisi mean field improved coupling we 
find 

— - = 4.22(2) — ► r A m = 0.625(3) (3.3) 

which is consistent with the value in eqn([3.2[), demonstrating that this coupling scheme is a rea- 
sonably good one. By contrast if we use a fit with the unadorned lattice bare coupling, g\{a), we 
find tqA^ = 0.541(3) which demonstrates that this is not a good coupling scheme. We can also 
modify the mean field coupling scheme by replacing the true value of the plaquette in eqn(2.2) 
with its perturbative expansion up to j-loops. We call this coupling scheme g 2 . . These Ij schemes 
will all have the same A parameter (since this depends on a 1-loop relation) however we find they 
work much less well than the / scheme. For example, the 1-loop improved coupling, I\, gives a fit 
leading to tqA^ = 0.448(2) - even worse than the bare lattice scheme! 

There is also a way to compare schemes directly, without needing an extra physical quantity 
like a/ro(a). This has the advantage that one can perform comparisons deeper into weak coupling. 
For a scheme s define the 3 -loop perturbative factor 

Now we expect for the SF scheme 

laA SF = {\ + C ^}F 3 SF [g 2 F (al)] (3.5) 
and for a lattice improved scheme 

oA/ = { 1 + c'a 2 } Fj [g] (a)] (3.6) 
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up to the various higher order collections. If we now replace the a 2 on the RHS of eqn(fT6|) by the 



expression for a in eqn(p.5|), and if we then take the ratio of the two equations, we obtain 

Asf _ lFf[g 2 SF (al)] {1 + ^} 



A, 



c = - 



I H\S 2 M] {l+c 2 ^{l + ^} 2 {^[^(a/)]F}" 



(3.7) 



We can now perform a fit for the constants c a ,c\ and c 2 over /3 ranges further and further into weak 
coupling, and see how rapidly Co approaches the known value of Asf/A/. In Fig. || we show a 
comparison for three schemes. Again we see that the Parisi scheme works well - and much better 
than the other schemes shown. 



R 




P 



Figure 4: Calculated values of Asp /As for the s — I, •, s = I3, o, and the s = L, x, lattice bare coupling 
schemes, all normalised to the known theoretical values. Horizontal errors indicate the range of /3 values 
used in each fit. 



4. Conclusions 

Taking advantage of the fact that large N lattice gauge theories have a well-defined weak 
coupling branch, we saw quite explicitly that 0(a 2 ) lattice spacing corrections are indeed important 
for transmuting the value of the bare lattice coupling into a value of the lattice spacing in 'physical' 
units [0]. 

We have also learned that the Parisi mean-field improvement scheme ^ for the bare coupling 
is in fact a reasonably good one. This we did by comparing it to the Schrodinger Functional scheme 
which we used as our benchmark. Obviously it will not be unique in this respect, and one could 
pursue this programme further. One cautionary remark: our benchmark SF coupling is defined in 
a finite volume, and one needs to understand the implications for this of the finite volume phase 
transitions at N = °° [|9|] that will lead to cross-overs at finite N. 



We can use fits of the form eqn(2T ) to extract values of Aj / y/o and hence A^/ y/a for all N. 
Doing so, in Fig. ||, we find that these values can be fitted with a modest 0(1 /N 2 ) correction 

A ^^0.503(2)(40) + °- 33 ^ 3 2 )(3) ; A>3 (4.1) 
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Figure 5: Calculated values of Aj^/ y^tJ versus 1 /N 2 with a linear extrapolation to N = °° shown. 

(We choose to exclude SU(2) from the fit, because of the difficulty in identifying a region where 
a weak coupling expansion is valid, but our fit does agree, when extrapolated to N = 2, with the 
value naively obtained there.) Here the first error is statistical and the second much larger error is 
expected to provide a bound on the systematic error from all sources. 
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